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Centrality dependence of charged particles multiplicity, transverse momentum spectra, integrated 
and differential elliptic flow, in y / s JVAr =200 GeV Au+Au collisions are analyzed using event by 
event ideal hydrodynamics. Monte-Carlo Glauber model of initial condition, constrained to re- 
produce experimental charged particle's multiplicity in 0-10% Au+Au collisions, reasonably well 
reproduces all the experimental observables, e.g. centrality dependence of charged particles multi- 
plicity, integrated and differential elliptic flow. Model predictions for higher flow harmonics, V3, «4 
however overestimate the experimental data, more in the peripheral collisions than in the central 
collisions. Correlation between initial (spatial) asymmetry measures and flow coefficients are also 
studied. With exception of the elliptic flow, for all the higher flow coefficients (?; n ,n=3-5), corre- 
lation is reduced with collision centrality. In peripheral collisions, higher flow coefficients are only 
weakly correlated to the asymmetry measures. Elliptic flow however, remains strongly correlated 
with initial eccentricity in all the collision centralities. 

PACS numbers: 47.75.+f, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 

It is expected that collisions between two nuclei at 
ultra-relativistic energies will lead to a phase transition 
from hadrons to the fundamental constituents, quarks 
and gluons, usually referred to as Quark-Gluon-Plasma 
(QGP). Experiments at the Relativistic Heavy Ion Col- 
lider (RHIC) at 75^=200 GeV Au+Au collisions [J 
00 and at the Large Hadron Collider (LHC) at 
7^=2.76 TeV Pb+Pb collisions 0001 had pro- 
vided compelling evidences for production of QGP. One 
of the experimental observables of QGP is the azimuthal 
distribution of produced particles. It is best studied by 
decomposing it in a Fourier series, 
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(1) 

4> is the azimuthal angle of the detected particle and ip 
is the plane of the symmetry of initial collision zone. 
In -yi^jY^OO GeV Au+Au collisions, second flow har- 
monic (1*2) , usually referred to as the elliptic flow, has 
been extensively studied experimentally as well as theo- 
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retically. Experimentally observed finite, non-zero v% is 
now regarded as definite proof of collective QCD mat- 
ter creation in Au+Au collisions. Qualitatively, elliptic 
flow is naturally explained in a hydrodynamical model, 
rescattering of secondaries generates pressure and drives 
the subsequent collective motion. In non-central colli- 
sions, the reaction zone is asymmetric (almond shaped), 
pressure gradient is large in one direction and small in 
the other. The asymmetric pressure gradient generates 
the elliptic flow. As the fluid evolve and expands, asym- 
metry in the reaction zone decreases and a stage arises 
when the reaction zone become symmetric and system no 
longer generates elliptic flow. Elliptic flow is early time 
phenomena. It is a sensitive probe to, (i) degree of ther- 
malisation, (ii) transport coefficient and (iii) equation of 
state of the early stage of the fluid. 

Ideal and viscous hydrodynamic models have been ex- 
tensively used to analyze the experimental data at RHIC 
and LHC energy collisions. Most of the analyses were 
performed with smooth initial matter distribution ob- 
tained from geometric overlap of density distributions of 
the colliding nuclei. For smooth matter distribution, the 
plane of symmetry of the collision zone coincides with the 
reaction plane (the plane containing the impact parame- 
ter and the beam axis). The odd Fourier coefficients are 
zero by symmetry. One of the important realization in re- 
cent years, is that the participating nucleons, rather than 
the reaction plane, determines the symmetry plane of the 
initial collision zone 0. The realization is the results of 
analysis of various experimental data, e.g. the two parti- 
cle correlation in A(j>- An plane [l(| [H| ■ The peculiar 
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structure in two particle correlations known as 'ridge' and 
'shoulder', observed both in STAR and PHENIX exper- 
iments have most compelling explanation provided the 
third flow harmonic, the triangular flow V3 develops in the 
collisions. Specifically, if initial condition is parameter- 
ized with quadrapole and triangular moments, response 
of the medium to these anisotropies is reflected in the 
two body correlation as ridge and shoulder [l3| , [l4| • Im- 
portance of the higher order flow harmonics in explaining 
the peculiar structures in two body correlation was also 
argued by Sorensen [l5[ . The ridge structure in pp colli- 
sions [lH [TtJ also has a natural explanation if odd har- 
monic flows develop. Recently, ALICE collaboration has 
observed odd harmonic flows in Pb+Pb collisions [l8[ . In 
most central collisions, the elliptic flow (i^) and triangu- 
lar flow (1)3) are of similar magnitude. In peripheral col- 
lisions however, elliptic flow dominates. More recently, 
PHENIX collaboration [Ig[ [20j] [2l| measured triangular 
flow in -^^^=200 GeV Au+Au collisions. 

In the present paper, in event-by-event hydrodynam- 
ics, with Monte-Carlo Glauber model initial energy den- 
sity distribution, we have simulated ^/s NN —200 GeV 
Au+Au collisions in 0-10% to 40-50% collision central- 
ities. Simulation results compare well with the exist- 
ing experimental data on charged particles multiplicity, 
transverse momentum spectra at pt < lGeV, integrated 
and differential elliptic flow. Higher flow coefficients how- 
ever are over predicted. We have also studied the cen- 
trality dependence of the correlation between the (inte- 
grated) flow coefficients with the initial spatial asymme- 
try measures. In 0-10%-40-50% collisions, elliptic flow 
remains strongly correlated with the initial eccentricity. 
Triangular flow is strongly correlated with initial trian- 
gularity only in very central collisions. The correlation is 
reduced significantly in peripheral collisions. Higher flow 
coefficients v n , n=4-5, even in central collisions, is only 
weakly correlated with initial asymmetry measures, e n , 
n=3-5 and the correlation is more reduced in peripheral 
collisions. 

The paper is organized as follows; in section UU the 
Monte-Carlo Glauber model for initial energy density for 
use in hydrodynamic simulations is briefly discussed. In 
section IIII1 hydrodynamic equations, initial conditions, 
equation of state used in the simulations are described. 
Results of the simulations are described in section IIVI 
Finally, all the results are summarized in section |Vl 



II. MONTE-CARLO GLAUBER MODEL OF 
INITIAL ENERGY DENSITY DISTRIBUTION 

In theoretical simulations of event-by-event hydrody- 
namics, one generally uses the Monte-Carlo Glauber 
model to obtain the initial energy density distribution in 
an event. Details of the Monte-Carlo Glauber model can 
be found in In a Monte-Carlo Glauber model, ac- 

cording to the density distribution of the colliding nuclei, 
two nucleons are randomly chosen. If the transverse sep- 



aration between the two nucleons is below ^2™, they 
are assumed to interact. Where <tnn is the nucleon- 
nucleon interaction cross section, taken here as 42 mb 
for Au-Au collisions at ^Jsnn = 200 GeV. Transverse 
position of the participating nucleons is then known in 
each event. The positions will fluctuate from event-to- 
event. If a particular event has N part participants, partic- 
ipants positions in the transverse plane can be labeled as, 
(xt,yi), {x 2 ,y2)--{x NpaTt ,y Npart ). Energy density distri- 
bution in the particular event can be obtained by as- 
suming that each participant deposit energy £0 m the 
transverse plane, 

AW* 

e(x,y) e Q ^2 6(x - Xi,y - yi) (2) 
i=i 

Fluid dynamical models require continuous density dis- 
tribution and discrete distribution as in EqJ5] cannot be 
evolved in a hydrodynamical model. To use in a hy- 
drodynamic model, the discrete density distribution has 
to be converted into a smooth energy-density distribu- 
tion. This can be done by smearing the discrete partic- 
ipating nucleon positions by some smoothing function, 
S(x -Xi,y- yt) -> g(x - x u y - y h &, ( 2 ~), d being pa- 
rameters of the smoothing function g. 

s{x,y)=e ^2 g(x-y,Xi,y-yi,Ci,C2--) (3) 
i=i 

One generally uses a Gaussian smoothing function. 
However, there can be other choices, e.g. in (23|, a 
Woods-Saxon distribution function was used for the 
smoothing. In the present simulations, we have used a 
Gaussian distribution 

(x-x i ) 2 + ( H - Hi ) 2 

9gauss{x - Xi,y - yi,a) oc e ^ , (4) 

of width cr=0.5 fm. 



III. HYDRODYNAMIC EQUATIONS, 
EQUATION OF STATE, INITIAL CONDITIONS 

With the Monte-Carlo Glauber model initial condi- 
tion for initial energy density, space-time evolution of the 
fluid, in each event, is obtained by solving the energy- 
momentum conservation equations, 

r"" = {e+p)u^u u -Pg^, (5) 
dfjT^ u = 0, (6) 

where e and p are the energy density and pressure re- 
spectively, u is the hydrodynamic 4-velocity. We have 
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FIG. 1: (color online) The black circles are the PHENIX data 
for the centrality dependence of charged particles multiplicity 
in y/s NN =200 GeV Au+Au collisions. The red symbols are 
the multiplicity in simulated events. The lines are to guide 
the eye. 



assumed ideal fluid formation and disregarded any dis- 
sipative effect. Assuming boost-invariance, hydrody- 
namic equations are solved in (r = y/t 2 — z 2 , x, y, rj s — 
^lnf^f) coordinate system, with the code AZHYDRO- 
KOLKATA [3. 



Hydrodynamics equations are closed with an equation 
of state (EoS) p — p{e). Currently, there is consensus 
that the confinement-deconfinement transition is a cross 
over. The cross over or the pseudo critical temperature 
for the quark-hadron transition is T c wl70 MeV (H-IH- 
In the present study, we use an equation of state where 
the Wuppertal-Budapest [1H H3 lattice simulations for 
the deconfined phase is smoothly joined at T = T c = 174 
MeV, with hadronic resonance gas EoS comprising of all 
the resonances below mass m res —2.5 GeV. Details of the 
EoS can be found in (29j . 

In addition to the initial energy density for which we 
use the Monte-Carlo Glauber model, solution of hydrody- 
namic equations requires to specify the thermalisation or 
the initial time and fluid velocity (v x (x, y), v y (x, y)) at 
the initial time. A freeze-out prescription is also needed 
to convert the information about fluid energy density and 
velocity to invariant particle distribution. We assume 
that the fluid is thermalized at Ti=0.6 fm and the ini- 
tial fluid velocity is zero, v x {x,y) — v y (x,y) — 0. The 
freeze-out temperature is fixed at Tp=130 MeV. We use 
Cooper-Frye formalism to obtain the invariant particle 
distribution of 7r~ from the freeze-out surface. Resonance 
production is included. Considering that pions constitute 
~ 20% of all the charged particles, tt~ invariant distri- 
bution is multiplied by the factor 2 x 1.2 to approximate 
the charged particle's invariant distribution. From the 
invariant distribution, harmonic flow coefficients are ob- 
tained as [301, 




FIG. 2: (color online) Charged particles transverse momen- 
tum distribution with and without multiplicity fluctuations. 



v n (y)e m ^ = 



dypTdpTd<f) 



fd<j)e mr t> ^ 

~ dN 

dyp T dp T 

dN 
dy 



(7) 
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In a boost-invariant version of hydrodynamics, flow co- 
efficients are rapidity independent. Present simulations 
are suitable only for central rapidity, y «0, where boost- 
invariance is most justified. Hereafter, we drop the ra- 
pidity dependence. ip n in Eqs l7l8l is the participant plane 
angle for the n-th flow harmonic. We characterise the 
asymmetry of the initial collision zone in terms of vari- 
ous moments of the eccentricity [13| , [14| , [3l| , 
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n = 2,3,4,5 (9) 



where x = rcos<p and y = rsin<f>. Eq|9] also determine 
the participant plane angle ip n . Asymmetry measures, £2 
and 63 are called eccentricity and triangularity. £4 and 
65 essentially measures the squareness and five-sidedness 
of the initial distribution. In the following, £4 will be 
called rectangularity. In the same vein, £5 will be called 
pentangularity. Fourth flow coefficient V4 is generally 
referred as hexadecpolar flow. In following, we refer it as 
the rectangular flow, which is more appropriate. v§ will 
be referred as the pentangular flow. 



IV. RESULTS 

A. Centrality dependence of charged particles 
multiplicity and pr spectra 

We have simulated 0-10%, 10-20%, 20-30%, 30-40% 
and 40-50% Au+Au collisions at ^^=200 GeV. In 
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FIG. 3: (color online) Black circles are PHOBOS mea- 
surements for the centrality dependence of elliptic flow in 
V^AW = 2fJ0 GeV Au+Au collisions. The simulation results 
for event averaged elliptic flow are shown as the red circles. 
The red squares, up triangles and down triangles are simula- 
tion results for triangular flow (1)3), rectangular flow (1)4) and 
pentangular flow (vs) respectively. 



FIG. 4: Event averaged flow coefficients (v n ) against the 
asymmetry measures, (e») for n=2-5. 



initial fluid velocity, freeze-out temperature etc. fit to 
charged particles pt spectra may be improved. 



each collision centrality, we have simulated iV el , e „ t =1000 
events. The constant £0 is fixed to reproduce experimen- 
tal charged particles multiplicity in 0-10% collision. It 
was then kept fixed for all the other collision central- 
ities. In FigfJJ simulated charged particles multiplici- 
ties are compared with the PHENIX data (32[. Once 
the model parameters are fixed to reproduce experimen- 
tal multiplicity in 0-10% collision, event-by-event simu- 
lations well reproduces the experimental multiplicity in 
other collision centralities. We do note that in peripheral 
collisions, simulated multiplicity overestimate the exper- 
imental multiplicity by ~10-15%. 

Even though charged particles multiplicities are well 
reproduced, the model simulations failed to reproduce 
charged particles pr spectra, in particular in the high pt 
region. In Fig|5J model simulations for charged particles 
p T spectra, in 0-10%, 10-20%, 20-30%, 30-40% and 40- 
50% collision centralities are compared with the PHENIX 
measurements 33]. Simulated spectra explains the ex- 
perimental data only up to pr ~1 GeV. In all the colli- 
sion centralities, at higher pr, model produces less parti- 
cles than in experiment. The results are interesting. It is 
well known that, compared to smooth hydrodynamics, in 
event- by- event hydrodynamics, pt spectrum is hardened 
[3C| . Still the hardening is not enough to produce requi- 
site number of particles at large px- It is also well known 
that pt spectra is hardened in viscous fluid. Better fit 
to charged particles pt spectra at large pt is expected if 
viscous rather than ideal fluid is formed in the collisions. 
We do note that in the present simulations, we have not 
made any conscientious attempt to fit the pt, spectra. 
Only the charged particles multiplicity in 0-10% collision 
was fitted. Varying other parameters e.g. initial time, 



B. Centrality dependence of flow coefficients 

1. Integrated flow coefficients 

Integrated flows are one of the important observables 
in heavy ion collisions. As discussed earlier, initial spa- 
tial symmetry is converted into momentum asymmetry, 
which is quantified in terms of different flow coefficients. 
For example, elliptic flow (t^) is response of an initial ec- 
centricity (£2) of the collision zone. Triangular flow (v^) 
is response of initial triangularity (£3) of the medium. 
Similarly, higher flow coefficients 1*4 and W5 are response 
of initial rectangularity (£4) and pentangularity (es)of the 
initial medium. In FigEl the black circles are PHOBOS 
measurements |34j | for centrality dependence of elliptic 
flow (112) ■ v% increases rapidly as the collisions become 
more and more peripheral. Present simulations for 1*2 in 
event-by-event hydrodynamics are shown as red circles. 
The simulations results agree well with the experimental 
data. In Fig|3l simulation results for (integrated) trian- 
gular flow (173), rectangular flow (174) and pentangular 
flow (1)5) are also shown. Triangular flow also increases 
as the collisions become more and more peripheral. How- 
ever, rate of increase is much slower than that for elliptic 
flow. V4 and v$ on the otherhand appears to be approxi- 
mately independent of the collision centrality. From cen- 
tral 0-10% to peripheral 40-50% collisions, they change 
by less than a few percent. 

In smooth hydrodynamics, elliptic flow in Au+Au 
collisions has been investigated in detail. Approxi- 
mately elliptic flow is proportional to initial eccentric- 
ity £2. Dependence of the event averaged flow coeffi- 
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: (a) 0-10% 



FIG. 5: (color online) In five panels, event-by-event simula- 
tions for flow coefficients v n , n=2-5 in Au+Au collisions are 
shown. The black, red, green and yellow lines are respectively 
for elliptic flow V2, triangular flow vz, rectangular flow 1*4 and 
pentangular flow vs . The black circles, red squares and green 
triangles are PHENIX measurements for elliptic, triangular 
and rectangular flow in Au+Au collisions at RHIC. 



cients ((«„), n=2-5) on the asymmetry measures ((e„)) in 
event-by-event hydrodynamics is shown in Fig01 The 
symbols, from left to right corresponds to 0-10%, 10- 
20%, 20-30%, 30-40% and 40-50% Au+Au collisions. As 
expected, asymmetry measures increases with collision 
centralities. The increase is most in £2, by a factor of 
~3.5 from 0-10% collision centrality to 40-50% central- 
ity. In other asymmetry measures, e n , n=3-5, the in- 
crease is more modest, factor of ~ 2-2.5 only. As it is 
in smoothed hydrodynamics, in event-by-event hydrody- 
namics also, elliptic flow increase, approximately linearly, 
with the initial eccentricity, (1*2) oc (€2)- Higher flow co- 
efficients, V3 also increase with initial triangularity, how- 
ever, the increase evidently is not linear. Still higher 
flow coefficients 1*4(1*5), approximately remains the same 
in all the collision centralities (as already shown in Figj3]), 
they appear to be independent of the asymmetry mea- 
sures, 64(65). Approximate centrality independence of 
higher flow coefficients, 1*4 and 1*5 in event-by-event hy- 
drodynamics indicate that unlike the elliptic or triangular 
flow, rectangular flow 1*4 or pentangular flow 1*5 may not 
be related to initial asymmetry measure of the collision 
zone. Later, we will discuss the issue in more detail. 
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FIG. 6: Correlation between elliptic flow (V2) and initial ec- 
centricity (62). Simulation results for V2 is plotted against 
the initial eccentricity for N event = 1000 events. For a perfect 
correlation V2 oc £2, all the points should he on a straight line. 



C. Differential flow coefficients 

Differential flow coefficients are very sensitive observ- 
ables and a model is well tested by comparing its pre- 
dictions against experimental differential flow data. In 
FigEl m nve panels (a)-(e), event-by-event hydrodynamic 
simulations for the differential flow coefficients, in 0-10%, 
10-20%, 20-30%, 30-40% and 40-50% Au+Au collisions 
are shown. In each panel, the black, red, green and yellow 
lines are the simulation results for elliptic flow, triangu- 
lar flow, rectangular flow and pentangular flow respec- 
tively. We have shown only the event averaged values. 
In each collision centralities, PHENIX measurements 
19] [20] 21] for the elliptic, triangular and rectangular 
flow are shown as the black circles, red squares and yel- 
low triangles. Simulations do reproduce the trend of the 
data, 1*2 > 1*3 > t*4 . Event-by-event hydrodynamic sim- 
ulations for the differential elliptic flow in Au+Au colli- 
sions agree well with the PHENIX data in all the collision 
centralities. In peripheral collisions, at pt > 2 GeV, ex- 
perimental data are marginally over predicted. We have 
simulated Au+Au collisions in the ideal fluid approxima- 
tion. If instead of ideal fluid, viscous fluid is produced, 
better agreement with data is expected. Indeed, explicit 
event-by-event hydrodynamic simulations 35] [36] 37] do 
indicate that the event averaged flow coefficients reduces 
with viscosity. In any case the agreement with data for 
elliptic flow measurements is much better for event-by- 
event ideal hydrodynamics compared to ideal hydrody- 
namic calculations with smooth CGC/Glauber model ini- 
tial conditions 38] 39] . 

Even though simulation results for elliptic flow reason- 
ably well agree with the PHENIX experiment, Simula- 
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FIG. 7: Same as in Fig[6] but for triangular flow (113) and 
initial triangularity (£3). 



FIG. 8: Same as in Fig[6] but for rectangular flow (114) and 
initial rectangularity (64). 



tion results for triangular (vs) and rectangular flow (V4,) 
appear to over predict the PHENIX data for the same. 
Interestingly, triangular flow is more over predicted than 
the rectangular flow. Also, the discrepancy between sim- 
ulations and experiment is more in peripheral collisions 
than in central collisions. For example, in 0-10% colli- 
sion, simulated V3, in the pr range 1-2 GeV, over predict 
the PHENIX data by ~ 30%. In 30-40% collision, the 
data are over predicted by ~60% or more. In 30-40% 
collisions, rectangular flow, in the pr range 1-2 GeV, is 
overpredicted by 5-10% only. Here again, better agree- 
ment with data is expected if instead of ideal fluid, vis- 
cous fluid is formed in Au+Au collisions. 

In FigJSl the yellow lines represents the pentangular 
flow 1)5. We did not find any experimental data for the 
pentangular flow. Following the trend of the simulation 
results for higher harmonic V3, U4, which are overpre- 
dicted in simulations, we do expect that the present simu- 
lation also over predict v$ . Present simulations then sug- 
gest that in experiments, in a peripheral 40-50% Au+Au 
collisions, in the px range 1-2 GeV, ~2-5% or less pen- 
tangular flow may be expected. 



D. Correlation between (integrated) flow 
coefficients and initial asymmetry measures 

Recently, in 40| [4l|. correlation between integrated 
flow coefficients (v n ) and initial asymmetry measures (e n ) 
of the collision zone was studied in event- by- event hydro- 
dynamics. It was shown that while elliptic flow remain 
strongly correlated with initial eccentricity, correlations 
between the higher flow coefficients v n and initial asym- 
metry measures e n , n=3,4,5, are much more weak. In 



40]|4lJ correlations between flow coefficients and asym- 
metry measures, in a single collision centrality, were stud- 
ied. How the correlations are affected, as a function of 
collision centralities were not studied. 

In Fig[Bl we have plotted the simulated elliptic flow 
(v 2 ) against the initial eccentricity (e 2 ) in N event =1000 
events. If V2 is perfectly correlated with £2, all the points 
should lie on a straight line. One observes that in central 
collisions, elliptic flow is strongly correlated with eccen- 
tricity. The correlation is marginally reduced in more 
peripheral collisions. One can conclude that in event- 
by-event hydrodynamics also, the correlation between el- 
liptic flow and initial eccentricity is strong, irrespective 
of the collision centrality. In FigfTJO results obtained 
for higher flow harmonics are shown. In central colli- 
sions, correlation between triangular flow (1)3) and ini- 
tial triangularity (£3) is strong, though degree of corre- 
lation appear to be less than that in elliptic flow. In 
more peripheral collisions however, correlation is signifi- 
cantly reduced. Correlation between higher flow harmon- 
ics, V4(vs) and initial asymmetry measure £4(65) even in 
central collisions is visibly much weaker than the corre- 
sponding correlation between elliptic flow and initial ec- 
centricity. The correlations deteriorate as the collisions 
become more and more peripheral. Indeed, from the scat- 
ter plot of V4 and v$ in peripheral collisions, it is difficult 
to claim that the flow coefficients are correlated with the 
asymmetry measures. 

In a quantitative measure was defined to quan- 
tify the correlation between flow coefficients and initial 
spatial asymmetry measure. A modified form is used 
here to quantify the correlation. For a perfect correla- 
tion, v n oc £„ and simulated flow coefficients will fall on 
a straight line. Dispersion of the flow coefficients around 
the best fitted straight line then gives a measure of the 



7 



(a) 0-10% 



11 - • • 

0.0 0.2 




0.01 ~, . 

ii>. 

0.00 ^W^^KLa^ 



40.0 0.2 



0.0 0.2 0.40.0 0.2 0.4 

0.02 



0.0 0.2 0.4 

e 5 



FIG. 9: Same as in Fig[6]but for pentangular flow (vs) and 
initial pentangularity (es). 



correlation. We thus define a correlation measure func- 
tion Cyneasurei 



measure 



y.\v 



random 



(e) - V st .line(e)] 



(10) 



Cmeasure essentially measures the dispersion of the 
simulated flow coefficients from the best fitted straight 
line, relative to completely random flow coefficients. It 
varies between and 1. If flow coefficients are perfectly 
correlated then v n oc e„ and C measure is identically unity. 
For completely random flow coefficients, C measure =0. To 
obtain an even ground for comparison of C measure for dif- 
ferent flow coefficients, the flow coefficients (v n ) and the 
asymmetry parameters (e„) are scaled to vary between 
and 1. In FigfTUl we have shown the correlation mea- 
sures for the flow coefficients as a function of collision 
centrality. The elliptic flow remain strongly correlated 
with initial eccentricity (c measur . e (t>2) ~ 0.95 — 0.99) in 
0-50% collisions. In central, 0-10%, 10-20% collisions, 
triangular flow (v^) is also strongly correlated with ini- 
tial triangularity (£3), (c measure (v2) ~ 0.95). Correlation 
is significantly reduced in more peripheral collisions and 
in 40-50% collisions, c measure (v2) ~ 0.75. In higher flow 
coefficients, correlation is even less in peripheral colli- 
sions. 

If departure of C measure from unity is interpreted as a 
measure of flow uncorrelated with the initial asymmetry 
measure, for elliptic flow V2, in 0-50% collisions, less than 
~ 5% of the flow is uncorrelated with initial eccentricity. 
In higher flow coefficients, v n , n=3-5, uncorrelated flow 
grow with collision centrality. For example, in rectangu- 
lar flow V4, uncorrelated flow grows from ~10% in 0-10% 
collisions to ^40% in 40-50% collision. 
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FIG. 10: (color online) collision centrality dependence of the 
correlation measure (see text) for V2, V3, V4 and V5 



As it was discussed previously, present analysis of 
Au+Au data also indicate that better description to the 
data will be obtained if viscous fluid rather than ideal 
fluid is formed in Au+Au collision. In |4(| viscous ef- 
fects on the correlation between elliptic flow and initial 
eccentricity and between triangular flow and initial tri- 
angularity were studied. It was shown that the correla- 
tions reduce significantly in viscous fluid. In more real- 
istic event-by-event hydrodynamic simulation of Au+Au 
collisions, with viscous fluid, correlation between higher 
flow coefficients and asymmetry measures then expected 
to reduce even more than obtained presently. 

Compartively low C measure or equivalently, large un- 
correlated higher flow harmonics is difficult to under- 
stand. In ideal hydrodynamics, final flow coefficients are 
related to the initial transverse energy density, or more 
appropriately on the pressure gradients only. We have 
even assumed zero initial velocity. Yet, the flow coeffi- 
cients, U3, V4 and V5 in peripheral collisions are largely 
unrelated to the initial asymmetry measures. Elliptic 
flow on the other hand is perfectly correlated with initial 
asymmetry measure. Source of the uncorrelated flows in 
higher harmonics can not be discerned presently. Possi- 
bly, other aspects of the initial density (higher moments 
or products of moments) are important in development 
of higher harmonics. 



V. SUMMARY AND CONCLUSIONS 

To summarize, in event- by- event hydrodynamics, we 
have studied the centrality dependence of charged par- 
ticles multiplicity, px spectra, and flow coefficients (in- 
tegrated and differential) in -^5^^=200 GeV Au+Au 
collisions. The initial energy density distributions are 
obtained from Monte-Carlo Glauber model. The Monte- 
Carlo Glauber model participant positions are smoothed 
with a Gaussian distribution of width er=0.5 fm and nor- 
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malized to reproduce experimental charged particles mul- 
tiplicity in 0-10% collision. We have simulated a large 
number of events, N event =1000 in each collision central- 
ity. Once the initial transverse energy density is fixed 
to reproduce multiplicity in 0-10% collision, the model 
reproduces the experimental multiplicity in other colli- 
sion centralities within reasonable accuracy. Experimen- 
tal charged particles transverse momentum spectra, how- 
ever are reproduced in the model, only in a limited pt 
range, pr <lGeV. At pt > 1 GeV, simulated spectra 
under predict the experiment. In the simulations, dissi- 
pative effects are neglected. Dissipative effect like (shear) 
viscosity, will enhances particle production, more at high 
Pt than at low pt- Better description to the pr spec- 
tra is expected if instead of ideal fluid, viscous fluid is 
produced in Au+Au collisions. We have also compared 
the model simulations for integrated and differential flow 
coefficients with experimental data. Experimental (inte- 
grated) elliptic flow in Au+Au collisions are correctly re- 
produced in simulations. The model also reasonably well 
reproduces the experimental differential elliptic flow in 
0-10%-40-50% collisions. In peripheral collisions, elliptic 
flow data however is overpredicted at high px- Higher 
flow coefficients and however are over predicted, 
more in peripheral than in central collisions. Here again, 
better description to the data is expected if instead of 
ideal fluid, viscous fluid is produced in Au+Au collisions. 



We have also studied correlation between (integrated) 
flow coefficients and initial asymmetry measures of the 
collision zone. In all the collision centralities (0-10% to 
40-50%) elliptic flow is strongly correlated with the ini- 
tial asymmetry measure, the eccentricity of the collision 
zone. The higher flow coefficients however show much 
less correlation with the corresponding asymmetry mea- 
sures. We have quantified the correlation and observe 
that with the exception for elliptic flow, which remain 
strongly correlated in all the collision centralities, for 
the higher flow coefficients, U3, U4 and V5, correlation 
reduces significantly in more peripheral collisions. It ap- 
pears that in higher flow coefficients, a significant part 
of the flow is unrelated to the initial asymmetry mea- 
sures. The reason for the flow unrelated to the initial 
asymmetry can not be discerned presently. One can only 
conclude that apart from the initial density distribution 
of collision zone, other aspects, e.g. higher moments or 
product of higher moments are also important for the 
development of higher harmonics. 
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